Dynamical Perturbation Theory 
and the Auxiliary Field Method 

B. Holdom* 

, Department of Physics, University of Toronto 

O's Toronto ON M5S1A7, Canada 

o\ ■ 

Qing Wang f 

China Center of Advanced Science and Technology (World Laboratory) 



(N 



Oh! 



P.O.Box 8730, Beijing 100080, China 
Institute of Modern Physics, Tsinghua University, Beijing 100084, China * 
Institute of Theoretical Physics, Academia Sinica, Beijing 100080, China 
Department of Physics, Tsinghua University, Beijing 100084, China 



■ Abstract 

ON 1 

By introducing auxiliary fields, which by their quantum numbers vanish in pertur- 



bation theory, we relate the dynamical perturbation theory of Pagels and Stokar and 
a successful gauged nonlocal constituent quark model to a U (N) gauge theory and to 



' QCD. This sheds light on the duality between quark models and resonance models. 

We then derive the effective action in the loop expansion without dropping cubic and 
quartic gluon couplings. 

1 Introduction 

Consider the full fermion propagator in the presence of a dynamically broken chiral symmetry 
in a gauge theory with no explicit fermion masses. 

S- 1 (p) = Z(p 2 )f,-Z(p 2 ) (1) 

S(p 2 ) vanishes to all orders in perturbation theory and may be considered to have the depen- 
dence e _1//fe9 on the gauge coupling g and some constant b. To account for this dynamical 
momentum-dependent mass function, Pagels and Stokar introduced "dynamical perturbation 
theory" (DPT) 0. We quote from their work. 
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To lowest order this approximation can be described as follows. Amplitudes 
that do not vanish to all orders of perturbation theory are given by their free- 
field values. For example, the amplitude Z(p 2 ) in (|I|) is given by Z(p 2 ) = 1. 
Amplitudes that vanish in all orders in perturbation theory such as A = e^ 1 ^ 9 , 
for example £(p 2 ), are retained, but we drop terms of order g n e~ l / bg , n > 0. The 
virtue of this approximation scheme is that amplitudes can be made to satisfy 

the chiral Ward identities consistently to lowest order in the parameters g 2 and 
A = e- 1 ^ 2 . 

One result of this work is a popular formula relating the (pseudo)Goldstone boson decay 
constant to the fermion mass function, which has been widely used in studies of QCD and 
other theories with chiral symmetry breaking. This formula is found to be in good agreement 
with a more sophisticated Bethe-Salpeter approach 0. 

Since the Ward identities are satisfied the idea of DPT at lowest order can be extended 
to a Lagrangian-based GNC (gauged nonlocal constituent) relativistic quark model |§ which 
preserves the chiral structure of QCD. The Lagrangian must be nonlocal to incorporate the 
momentum dependent mass. The GNC model reproduced the Pagels-Stokar decay constant 
formula, and when expanded to order p A in the energy expansion was found to be surprisingly 
successful at describing low energy chiral dynamics. The momentum dependence of the quark 
mass function provides a natural ultraviolet cutoff on the quark loop integrals, thus leading 
to the correct Wess-Zumino terms without regularization subtleties [|J. The incorporation 
of this momentum dependence represents an improvement over constant mass quark models 
of the Nambu-Jona-Lasino type. But the question remains, why does a free quark model 
represent low energy QCD dynamics as well as it does? For example at the level of the quark 
propagator, how can the momentum dependence of Z(p 2 ) in ([J) be ignored? 

By using the auxiliary field method we shall provide a systematic derivation of DPT and 
the GNC quark model from a gauge theory. The connection with DPT is made by introducing 
auxiliary fields only for those quark bilinears which define purely nonperturbative amplitudes. 
The auxiliary fields introduced are thus chirality changing, and they may also carry color or 
fermion number. The effective theory at lowest order in a loop expansion will correspond to 
the lowest order in DPT, and the origin of the GNC quark model is made clear. 

We shall also find a set of ladder SD equations which are all homogeneous equations in 
the auxiliary fields. It may be noted that these SD equations treat all the possible quark 
mass functions in the various channels, including for the usual qq mass, on an equal 

footing. In the SD equation for we will automatically have Z(p) = 1, which is usually 
obtained in studies using the ladder SD equation via a special choice of gauge. In our 
approach Z(p) = 1 follows at lowest order in the loop expansion. 

We shall deal with the problem of gauge boson self-interactions by studying a strongly 
interacting U(N C ) (rather than SU(N C )) gauge theory with Nf flavors of quarks. In U(N C ) 
there are two coupling constants, g and g\ for SU(N C ) and U(l) respectively. Both of 
these interactions are attractive in the usual color- and flavor-singlet qq channel. For chiral 
symmetry breaking to occur in this channel there is some constraint on some combination of 
g and g\. We shall take advantage of the fact that we can hold the scale of chiral symmetry 
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breaking fixed while g is free to vary. In this way the theory can smoothly interpolate between 
the QCD limit (large g and small gi) and the U(l) limit (small g and large gi). For some 
subset of this space of theories it makes sense to treat g perturbatively on the momentum 
scales which characterize chiral symmetry breaking. For some scale sufficiently below the 
chiral symmetry breaking scale g will grow strong and cause confinement. But the effects 
due to the nonabelian nature of SU(N C ), namely trilinear and quartic gauge couplings, may 
be treated perturbatively as far as the study of chiral symmetry breaking is concerned. 

Our approach will also shed some light on the concept of duality in the description of 
low energy chiral dynamics. This is related to the observation that the low energy chiral 
Lagrangian to order p 4 of QCD is surprisingly well described by constituent quark models, 
such as the GNC quark model. The quark degrees of freedom are somehow offering a dual 
description to the effects of the lowest lying vector and axial vector resonances, which are 
usually thought to dominate low energy chiral dynamics. We shed light on this by showing 
that QCD may be rewritten as a path integral over various auxiliary field degrees of freedom 
which do not include vector or axial vector degrees of freedom. Instead we will have scalar 
and tensor fields, some of which have color and/or diquark quantum numbers. The lowest 
order in a loop expansion involving all massive degrees of freedom corresponds to setting all 
the non-Goldstone-boson degrees of freedom in the various auxiliary fields to their vacuum 
values. The fact that we are not required to fix vector or axial-vector degrees of freedom 
then suggests how it is that a constituent quark model can reproduce the effects of vector 
or axial-vector degrees of freedom. 

Our basic result is that the perturbative expansion of a nonabelian gauge theory may be 
reorganized in such a way that the lowest order in the new loop expansion coincides with 
a previously studied and successful quark model of low energy chiral dynamics. Instead 
of being put in by hand, the quark mass function emerges from the dynamics. Also the 
quark model now exists in a framework of a loop expansion in which corrections, at least in 
principle, may be systematically worked out. Such a process was implied but not specified 
in the original proposal for DTP. As one example the phenomenological effects of various 
scalar and diquark degrees of freedom could be explored in the context of a quark model, 
without fear of double counting degrees of freedom. 

In the next two sections we shall derive generating functionals without and with, respec- 
tively, the effects of cubic and quartic gauge couplings. The latter case is significantly more 
complicated, and it will be carried out by introducing two sets of auxiliary fields. In the 
next section we also show how the (pseudo)Goldstone degrees of freedom may be retained 
for a description of low energy physics. 

2 Low energy theory 

Due to our study of U(N C ) rather than SU(N C ), we will first ignore effects from the cubic 
and quartic gauge couplings from the beginning. Then the gauge fields may be integrated 
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out to give the following generating functional. 



iW[J] 



2>02>0expz| I d 4 xi\){i(j) + J)ip 



2„2 



* 9 



2 

,'2„2 



A 



(2) 



G fiu (x,y) and Gfy(x,y) are the tree level SU(N C ) and 2/(1) gauge field propagators respec- 
tively. J represents the external sources which are nontrivial matrices in the SU(Nf) flavor 
space. 

J(x) = j) (x) + (x) 75 - s(x) + ip(x)7 5 (3) 
The generating functional may be rearranged as follows. 



iW[J] 



■2 2 

VipVipexpi^J d A xijj{i(j) + J)ip + J d A xd 4 yG '^(x, y) 
+ [?L,Ja;)(^) Q , 7 ^l i ,(x)][^ LiQ ,( 2/ )(^)^ 7 >i i/3 ,( 2/ )] 
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+ Wn, a (^UY^rAx)} W R , a < (y) (f W7>j^ (y)] 



d 4 xd 4 yG^(x,y) 



+ W L , a (?)-f1>U (*)] (y)7"^ (?/)] + fea (*) VVka (*)] 0** O/WL' (*/)] 



VipVipexpi y I d xip(i$ + J)ip 
1 



d 4 xd 4 y 



(4) 



where 

= ^^(^.y)(y)^(V i ) (rp )(Y)^(7' i )^ + ^! t 1 j5 Q/3 (^) (Tp )W(7 M ) (7V (5) 

We may follow the auxiliary field method || and insert a Gaussian integral to cancel the 
four-fermion terms. Note that our choice of auxiliary fields differs from previous treatments. 



constant = / Vxr^Xl^ KrV KjJD K R exipi J d 4 xd 4 y 
x {Tr [{xr ~ K^ R ?p L ) K- 1 ( X l ~ K^ R 



— Tr 
2 



(k l - K^^j) K- 1 (k l - Ki/j l i/j l ) + (L - R)] | 



(6) 



The various auxiliary fields are bilocal, xl{%,u), etc., and the traces imply contraction of 
indices similar to We will not need the explicit form of K~ l . The resulting action takes 
the form 

e iW[J] = J V X RVx L VK L VK R VK L VK R V^V^expi I^J d 4 x^(i^ + J)ip 
+ J d 4 xd 4 y[Tr (xrK- 1 Xl) ~ ^Tr {k^kl) - ^Tr {k r K' 

-^LXR^R - *pRXLlpL ~ IpL^LlpL ~ 1pL K Llp L ~ IpR^RlpR ~ iPrKr^r } (7) 

The fermion fields could now be integrated out. The original path integral over the 
quark and gluon fields has been traded in for a path integral over various bosonic degrees 
of freedom. We first consider the lowest order in the loop expansion with respect to all 
auxiliary fields. The stationary condition for the tree level action with J = gives a set of 
SD equations, homogeneous in the auxiliary fields. The auxiliary fields may be decomposed 
into various irreducible representations of the SU(N C ) x SU(Nf) color and flavor symmetries. 
If we define 

X = Xl + Xr (8) 

then the SD equation for the flavor- and color-singlet degree of freedom, x( x jV)%ij ~^ 
SapSpaSij^x - y), is 



E(x - y) 



2 m-r 



G^(x-y)+g 2 1 G^{x-y) 



xY [P L S(x - y)P L + P R S(x - y)P R ]Y, 



(9) 



with iS (x — y) = i$5(x — y) — H(x — y). (10) 

In the presence of an ultraviolet cutoff this equation will determine a mass function 

All other auxiliary fields are in less attractive channels and we may assume that they take 

on the vanishing solutions. 

We note that except for the (pseudo)Goldstone bosons (PGBs) associated with chiral 
symmetry breaking, all other auxiliary field degrees of freedom are expected to have masses 
of order the chiral symmetry breaking scale or higher. This leads us to consider the lowest 
order in a loop expansion with respect to all massive, non-Goldstone degrees of freedom. 
That is, we freeze all massive auxiliary fields to their vacuum (J = 0) values as determined 
by the SD equations above, but we wish to retain the light PGB degrees of freedom as 
dynamical fields in a low energy theory. Note that nonzero sources J are needed to define 
the generating functional of the low energy theory. Nonzero sources are consistent with the 
freezing of the massive fields at their J = vacuum values as long as the sources only vary 
on distance scales large compared with inverse masses of the massive degrees of freedom. We 
will say more about the effects of the massive degrees of freedom at the end of this section. 

We constrain the massive degrees of freedom as follows. 



ti c [\ a x{x, y)\ 
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tr[<vx(a:,y)] 



a ( x > y)ik a (yi x )kj + n 0> 2/L n (^ 



;r 







a = tr( X )/4, n 



V)' 
tr(i7sx)/ 4 



(12) 
(13) 
(14) 



tr (tr c ) denotes a trace over spinor (color) indices only. The first constraint keeps x{ x i y) 
diagonal in color space, and we henceforth ignore color indices. The second constraint sets 
the tensor degrees of freedom to zero. The third constraint, where i, j are flavor indices, 
allows for PGB degrees of freedom in flavor space while constraining the scalar mode. 

With these constraints on the fields we may now go back and consider the K~ x terms in 
(0). These terms are quadratic in a and 7r and by chiral symmetry they must take the form 



d 4 xd 4 y [a ij (x,y)a ji (y,x) + Ely (x, y)Hji{y, x)\ f{x - y) 



(15) 



The constraint (|13| ) then implies that these terms are independent of the PGB fields. The 
PGB degrees of freedom are then described completely by the Trln terms, and the generating 
functional for slowly varying sources is 



JW[J] 



U{x) 



J VUexp [-iTrln (i@ + J + Xit,j( x , y))) 



2iir(x)"f 5 1 fir 



(16) 
(17) 



X%,j{ x iV) must be chosen to satisfy the constraints above and to transform correctly 
under chiral transformations of 7r(x) and J(x) = f [x] + ^ (x)j5 — s(x) + ip(x)^. A choice 
for x-k,j{ x iU) corresponds to a choice of representation of the ir{x) fields. A minimal choice 
is the following ||. 



Xn,j(x,y) = S(x - y)[£(x)X(x,y)£(y)] 



X(x, y) = Pexp 



r M (z)c^ 



r, 



U(x) 
1 



(18) 
(19) 

(20) 

(21) 



The origin of the dependence on s(x) and p(x) is explained in 0. 

We have arrived at the GNC quark model, which when expanded in powers of derivatives 
and quark mass, successfully models the low energy chiral Lagrangian of QCD. But our 
derivation raises a further issue. When we froze all the massive auxiliary field degrees of 
freedom, we not only removed the effect of these fields in loops, but we also removed their 
tree-level effects. In the context of a loop expansion, at lowest order we should keep the 
effects in the low energy theory of the tree-level exchange of massive fields. In particular 
the original GNC model should be supplemented by the effects of scalar exchange. This 
would supplement the quark loop effects, which by our discussion of duality, are related to 
the effects of vector and axial-vector resonances. 
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Figure 1: Effective interactions caused by integrating out the gauge fields. The blob with 2n 
external legs stands for the full gauge boson n-point function without fermion loops. The 
small empty block stands for the fermion field. 

3 Effective Action 

In this section we shall present a derivation of the effective action and the loop expansion 
without dropping the cubic and quartic gluon couplings. The auxiliary fields we introduce 
will again correspond to nonperturbative amplitudes only. After integrating out the gauge 
fields, the generating functional becomes 



Z[J] 



,iW[J\ 



d 4 X^ + J) + E / ■ ■ ■ d^ n -K\^X{xu ■■■,Xn) 

n=2 J n - 

x[C(^i)C(^)] ■ • • 0:: n M4TM] (22) 

See Fig. 1. The Kalfa-^pi (x%, ■ ■ ■ ,x n ) are defined as follows 



T^cripl---cr n p n / 



+Sn2glG^ 2 (Xl , X 2 )8 aiPl {Y 1 )<r ipi <W 2 (l 



n 



Ian 



ICT2P2 



(23) 



G^"!^ i (xi, ■ ■ ■ , x n ) is the full gluon n-point function without fermion loops included and 
Gfy(x,y) is the U(l) gauge field propagator without fermion loops. 
We may rewrite the sum over n in (E21) as follows. 



oo „ 

/ d 4 xi ■ • ■ cf 

n=2 

x {xxWC Ml KT W1>r M\ ■ ■ ■ RT M4ZT- 1 (xn-i)] 

oo „ l-\2ns~in 

Jd Xx-'-d X 2n . ,j K-cah-aanlhn ^1' ' ^ 



n=l ' 
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Figure 2: The chirality (denoted by L and R) changing part is explicitly separated out, such 
that diagram la (from Fig. 1) is separated into diagram 2a (the chirality changing part) and 
diagram 2c, and diagram lc is similarly separated into diagrams 2b and 2e. 



X WrX {KlWlt: (**,)] K (*l)] • " " 

• • • w^::T\^ n -^^^ 

oo „ i-\n 

+ £ / d* Xl ■ ■ ■ Jx^KZgZfcixu ■■■,x n ) 

n=2 J n - 
XC^OC ^)\WZ\^T M ■ ■ ■ W:I n {XnWl Pn {Xn)\ (24) 

where K is defined by this expression and C™ = m \^L m y - See Fig. 2. Inserting this back 
into the generating functional and introducing auxiliary fields gives the following. 

Z\J\ = e^ J ] 

= j VS R VS L VFVFVBVBV^5 (sgjfo y) - K, a {xW^{y)) 

xS (f$»(x, y) - ^{x)^{y))5 (b^^ ( Xl ,x 2 , x 3 ) - ^ (*i)#T (* 2 fiC (* 3 )) 
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x5[B^™°*{x l ,X2,x 3 ) ~ CWCWCW ) ex PM ld 4 x^+J)^ 



+ 



ja xi---a x 2n A QiA ... a2n/32n ^i,---,x 2n j 



C ,4 2n-l°"2n-l«2n-2P2n-2 
i,a2n-l/32n-2 



f \ C i 2nO-2n«2n-lP2n-l / \ 

^2n-l, ;r 2n-2j^ J R )a2n /32n-i \ X 2n, X 2 n-l) 

*) 2 " 



oo „ 

- E y • • • ^ 4 ^^^^i::: 2 2 X(^> • • • , 



-f7r«2nO"2n*lO"l / s -F7«2n-2f2n-2*2n- 1 C«2n- 1 / \ 

X - t a 2n a 1 [X2n,Xl) • • " a2n _ 2 a 2 „-i (, x 2n-2 , ^2n-lJ 



v p«lPl«2P2/ N 



77«2n-lP2n-l«2nP2n / \ 

^/3 2 n-i/3 2 n l x 2n-l, x 2nj 
\2n+l 



«10"1« 2 (T2 130-3 



aia 2 «3 v X l> X 2, x 3j-T Q 4Q5 l x 4, x 5j • • •- fV Q2„a2n+l l x 2n, x 2n+lj 



r>llPl*2P2*3P3 
A - tJ ai02Q!3 



(xi,x 2 ,x 3 )i4gf t ' ,B (x4,X5) • • -^r^r iP2,l+1 (^,a; 2n+1 )} (25) 



See Fig. 3. We can remove the delta functions by instead including the following terms in 
the exponential and integrating over the additional auxiliary fields, xl, Xr-> U, U, V, and V. 



i I d 4 xd 4 y 



+ [Sl^{x, y) - i>R, a (x)i>L, a (y) ) xl^x, y) + ( ^T( x > y) - i>*{x)tfp{y) ) u^ p (x, y ) 
+u% p (x, y) [ KT^ v) - #r(* Wv) 



+i J d 4 Xid 4 x 2 d 4 x 3 



B 



+V. 



1CT1«20'2«30"3 / 



,1CT1«20'2«30"3 

il<Tli2<T2«3(T3 / 



(x 1 ,x 2 ,x 3 ) - ^ a 1 \x 1 )^ 2 {x 2 W^\x z ) )V C 



(x 1 ,x 2 ,x 3 ) 



See Fig. 4. We may write the result as follows. 

e lW[J] = JvS R VS L Vx R Vx L VFVFVUVUVmBVVVVexpi{D[x R ,XL,U,U,V,V,J} 

i S R,ap( x , y)XR™p(x, y) + SZ^x, y)x l ^p(x, y)} + T°p P (x, y)U l °l P {x, y) 



d 4 xd 4 y 



+U l lf) P {Xi y)K7 P ( x ' V) + d 4 x 1 d 4 X2d 4 x 3 [B^l^ 3 2t3a3 (x 1: x 2 , x 3 )V l a \ 



<71«20"2i30-3 / \ 

Ql a 2 a3 l x i, x 2,£3j 



I 12^213^3/ \ R il(Tli2CT2i3CT3^„ „ „ \1 

+ V {Xl,X 2 ,X 3 )tf aia2as {Xi,X2,X 3 )\ 



00 



+ E / ^ 4xi ' ' ' ^ 4;r 2r. 
n=l J 



102C«3 
\2n/^(n 



\V W 2n r r O-ipi-CT2nP2n / \ 
(2n)' a lft- a 2n/3 2 nl X l' " ' ' X 2«J 
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Figure 4: These new terms are caused by realizing the constraints among the fermion fields 
and the composite auxiliary fields. 



X(5 L,ai/3 2 „ \ X ^-> X ^njO Ra2f5l [X 2 ,X 1 ) 

Cf»2n-10"2n-l*2n-2P2n-2 / „ \C'2 

■ ■ ■ ^L,a 2n -ift>n-2 l X 2n- 1 , X 2n -2 ) Ci R 

\2n 



«2nC2n*2n-lP2n-l 
"2n/32n-l 



^2n, ^2n-l) 



- e / ^ • • • ^-^<x::rjr>i, • • • , **o 

-p«2nO"2n«10"l / s «2n-2 02n-2 *2»- 1 "2m- 1 

Xi< a2»«l V^ 2 ™' ^1/ ' " t «2n-2Q!2n-l 



v pilPli2P2 ( ' \ T?i2n-lp2n-li2nP2n( \ 



/3ift> 
oo „ 

n=l 



d xi---d x 2n +i 



(i) 2n+1 



(2n+ 1)! 

(Xi,X2,X3)F, 



^2n-25 ^2n-l) 
C2n) 

^CTlpi--CT2„+lP2n+l / \ 



/ \ -7=r*2nO"2n«2n + 10"2n + l / \ 

(x 4 , x 5 ) • • • F a2na2n+1 (x 2n , x 2n+1 ) 
xB%^"( Xl ,x 2 ,x 3 )F^»(x A ,xd ■ ■ •i J S^Sr ipa - +1 (x a „,x an+1 )} 

The path integration over the fermion fields has been isolated. 



expi£>[xi?, Xl, f 7 , Ef, V", V, J] = / VipVipexpi 



d 4 xd 4 y 



^L,cc{ x Wii}{y)x l R,ap( x ^ v) 



— / d 4 x\d 4 x 2 d 4 xz 



ao"3 



10-1120-2*30-3 , 
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+v 



l<Ti«2<T2«3<T3 



x 1 ,x 2 ,x 3 )^(x 1 )^(x 2 )^p(x 3 ) 



(27) 



We may use the loop expansion to complete the path integral over the fields xr-> Xl,U, 
U, V and V. 



iW[J] 



VS R VS L VFVFVBVB exp i\ T[ X r, c , X l,c, U c , U c , V c , V c , Sr, S l , F, F, B, B, J] 



(28) 



+^i[Xr,c Xl,c, U c , U c , V c , V c , J] 



where iT[xr, Xl, U, U, V, V, Sr, Sl, F, F, B, B, J] is the argument of the exponential in ( p6| 
and 



Tl [XR,c, XL,c, U c, U c , V c , V c , J] 



(29) 



-Trln 

2 



d 2 m<s>,j] 



all 1PI Feynman vacuum diagrams 



with inverse of propagator — : - 



and n-point vertex — i—rr- 



d n T[(f), $, J] 



xfa u( Pa n f3 n \ x n,X n ) 



n = 3,4 



(30) 



is the set of fields Xr,Xl,U ,U ', V ,V , and $ is the set of fields Sr,Sl,F,F,B,B^ Ti[<p c , J] is 
independent of $ since the T[<f>, $, J] in (|3~0|) can be replaced by D[(f), J] (since T[</>, $, J] — 
D[(f>, J] is linear in <fr). The classical <p c fields satisfy the equations 



d{T[ ( Pc^,J}+Ti[<Pc,J}} 

d<t> c a ,ap(x,y) 



(31) 



which determines <p c as a functional of $. We may again use a loop expansion to integrate 
out the field $. 

^[j] = r s [$ c ,j] + r sl [$ c ,j] (32) 



r s [$ c , J] 
r s i[$ c ,J] 



r[</> c ,$ c , J] + r a [0 c , J] 



-Trln 



a 2 r s [$,j] 



(33) 

+ all 1PI Feynman vacuum diagrams 



with inverse of propagator 
and n-point vertex 



d 2 r s [$,j} 



9 n r s [$, j] 



n = 3,4,.... 



(34) 



1 V", V, B, B should appear with three indices instead of two. 
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The $ c fields satisfy the equations 

d{T s [<t> c ,J] + T sl [^ c ,J]} 
V) 

which with the help of (^) and (j31~D become 

9{r[0 c ,$ e ,j] + r gl [$ c ,j]} 



(35) 



(36) 



Note how the results in (|31|) and (|3^) have a similar form. 

The fields U,U,V,V, F,F,B,B carry nontrivial fermion number and color, while fields 
Xl,Xr, Sl,Sr and external sources do not. Because of our dynamical arguments associated 
with the most attractive channel we may take the nonvanishing classical fields, for arbitrary 
J, to be colorless and fermion-number conserving. 

$ c = for $ C ?S L ,S R (37) 

<Pc = for <f) c ^ xl,Xr (38) 

S fi%^v) = ^Sf p {x,y) (39) 

Xf'a^v) = ^Xa JP (x,y) (40) 



L > 



The result is that only the first two diagrams in Fig. 3 and and the first four diagrams in 
Fig. 4 survive. Then the generating functional is 

W[J] = -iTrln(^+ J- PlXlPl-PrXrPr) 

+W S [S R , S L ) + TxIxr, Xl, J] + T sl [S R , S L , J] (41) 



with 

W S [S R ,S L ] = EJ* x i'--*xi n ^^KZ£;%£(x 1 ,-.. } x 9n 



n=l 



(2n) 

xsr ii2np2n (xi,x 2n )s i r iipi (x 2 ,x 1 ) ■ ■ ■ 

Qi2n-lC r 2n-li2n-2P2n-2 ( \ Q^n^nhn- 1 P2n- 1 / \ / ^ r)\ 

• • • \X2n-l, X2n-2jO R \X^m x 2n-l) 



If we apply ([H|) we find that 



[JW7(y, x)Pr\^ + y) + ^ d l ll Sf XL ^ ] = ° ( 43 ) 

7V C Oxr {x,y) 



-iS^(y, x) — [(0+ J — P lX lPl - PrXrPr)- 1 ) (v, x) (44) 
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On the other hand (p6|) becomes 
XR JP (x,y) = 



dw s [s R ,s L \ ar 8l [s R ,s L ,j] 



N c ldST p {x,y) dSr p (x,y) 



(45) 



By using the last two equations, (fy]) may be written as 

W[J] = tTi\nS + W s [S R ,S L }+r 1 [ XR ,XL,J}+r sl [S R ,S L ,J} 
— / d 4 xd 4: y 



nicrjpf ,9W S [S R , S L ] niujp, ,dW s [S R , S L ] 



dS^ p (x,y) 



dSr%x,y) 



nicrjp/ \9T s i[Sr, Sl, J] niajp, sdT s l[S R , Sl, J] 

-S Rc P {x,y)—- I — / Sri p (x,y) 



dS%/(x,y) 



dSZl P (x,y) 



(46) 



We may also rewrite (|44|) as 



N, 



f dW s [S R ,S L ] dT sl [S R ,S L ,J] \ 

•Oil or. 7 \ 1 \ K-R 



dS R (x,y) dS R (x,y) J 

+Pl r dw s [s R ,s L ] + dv al [s R ,s L ,J] \ ^ 

L \ dS L (x,y) dS L (x,y) J 1 



IV] p 



(47) 



These are our final full results, with ( f47l) playing the role of a generalized SD equation when 
used along with ([43]). 

If we consider results at lowest order in the loop expansion then 



W[J}=iTi\nS + W s {S R ,S L } 



d 4 xd 4 y[S l R JP (x, y) dW °l S *> Sl ) + St jp (x, y fW'&R'SL] 



tS- 1 - (i$ + J) 



1(7] p 



[x,y) 



PnS^(y } x)Pn + ST p (x } y)=Q 



dW s [S R ,S L ] dW a [S R ,S L ] 

r R , r ± R + f L or7 / —-TL 



dS R (x,y) 



dS L (x,y) 



to -]p 



(48) 

(49) 
(50) 



We notice that the dependence on xl and xr nas been removed. We could instead express 
all results in terms of xl and xr by using (f44"l). 

Alternatively we note the result ( pEo] ) which now reads 



Xr JP (x, y) 



1 dW s [S R ,S L ) 



(51) 



N c dS^ p (x,y) 

L 

This defines the relation between S and x fields and suggests a Legendre transform 
W x \xr, xl] = W S [S R , S L ] +N c f d*xd 4 y[S% jp (x, y)x i R jp (x, y) + Sr P (x, y)xT P (x, y)} (52) 
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N c SZ jp (x,y) (53) 



such that satisfies 

dW x [ X R,XL\ 
dx l R P {.x } y) r. 

L 

We may thus write our generating functional ([52]) as 

W[J] = -zTrln(^ + J-PlXlPl-PrXrPr) + W x [ X r,Xl} (54) 

To make contact with the last section, if we keep only the lowest order term in g in 
W s [Sr, Sl] then W s [Sr, Sj] is bilinear in the fields. As we described in the introduction, 
the neglect of these higher order terms could be justified in the U(N C ), but not the SU(N C ) 
theory. From ([42]), 

N r a 2 

W S [S R ,S L ] = J d 4 xd 4 y[^G a ^(x,y) + glG^(x,y)]Tr[YSl(x,y)YSl(x,y)} 

= -^/^[«] (55) 
In this compact notation x r — fcST, which with the help of (j50|) , and fl5"I|) can be 

L R 

written as 

X n(x,y) = - l l^G™(x,y) + g 2 G$(x,y)}rP LR (i@ + J-PLXLPL (56) 

This is the SD equation written in terms of the x fields. Finally we have 

at r 

W x [ X r,Xl] = -^Jd^yTrixLlC-'xR} (57) 



With this last result and (54) we essentially recover the form of the generating functional of 
the previous section, with a minor difference being that gluon loop effects are now included 
in the gluon propagator. 
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